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Abstract 

A study of zero- dimensional theories, based on exact results, is presented. 
First, relying on a simple diagrammatic representation of the theory, equa- 
tions involving the generating function of all connected Green's functions are 
constructed. Second, exact solutions of these equations are obtained for sev- 
eral theories. Finally, renormalization is carried out. Based on the anticipated 
knowledge of the exact solutions the full dependence on the renormalized cou- 
pling constant is studied. 
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1 Introduction 



In this paper we study several aspects of zero-dimensional quantum field theory. Such 
theories may serve as a model (the static ultra-local limit) of more realistic quantum 
field theories, and as a useful didactic object in their own right, since zero-dimensional 
theories, for which the path integral is actually a simple integral, allow for many explicit 
and exact solutions that cannot be obtained in higher dimensions. As recent examples, 
we may quote 't Hooft |L| and Bender et al. 0. Questions of particular interest here 
are the behavior of theories in high orders of perturbation theory (either many loops, or 
large number of external legs), and of the relation between the diagrammatic perturbation 
expansion and the full solution. The layout of this paper is as follows. We start by 
a diagrammatic (re) derivation of equations that govern the set of all connected Green's 
functions of the theory. We show how for a general scalar theory with arbitrary interactions 
the Green's functions may be obtained order by order. We point out how the Schwinger- 
Dyson equation, although derivable from purely diagrammatic arguments, in fact describes 
a much larger class of solutions. Next, we discuss the representation of these solutions as 
path integrals over contours in the complex </>-plane. Exact solutions for theories with 
interactions up to <£> 4 are obtained including explicitly perturbative and non-perturbative 
contributions. We show how a classification of the allowed contours in the complex <p- 
plane can immediately determine properties of the non-perturbative character of these 
theories. Renormalization, which for these theories is equivalent to imposing restrictions 
to diagrams, is also studied. The wave function renormalizations are fully determined and 
their dependence on the renormalized coupling constant of the theory is presented and 
discussed. 
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2 Basic equations 



In this section we derive equations for an arbitrary zero-dimensional field theory. The 
derivation is based entirely upon the diagrammatic representation of the theory. A theory 
is diagrammatically defined by a sequence of vertices that are weighted by the 'coupling' 
constants taken for convenience as — for the k-th vertex. In fact, the two-point coupling 
A2 = m 2 = fi can be eliminated by the introduction of the propagator which means that 
every line of a diagram accounts for a factor Moreover a loop in a graph is counted 
by an additional parameter, h. A solution of a zero-dimensional theory is determined by 
a sequence of objects C n , n = 0,1,2,... that represent the connected, n-point Green's 
functions, i.e. the sum of all connected diagrams with n external lines. One can define the 
generating function of the Green's functions as 



00 

x 



= E ■ (!) 



n=0 



We want to write down an equation for 0, and in order to do so, we represent it with a 
diagram: 

Derivatives of with respect to x are represented by extra lines, 

m = — o— > </>"(*) = 

and so on. 



2.1 The Schwinger-Dyson equation 

Let us consider a theory with a k-th vertex. In order to write an equation for we start 
with the bare vertex and we attach k — 1 blobs 

The factor l/(k — 1)! is due to the k — 1 identical blobs. Considering a one loop attachment 
we similarly have 



2! (k - 3)! ' 
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The factor l/(k — 2)! is again due to the k — 2 identical blobs, whereas the 1/2! is due to 
the symmetry factor of blob with two lines. Following this reasoning we can proceed with 
higher powers of H. For instance at the two-loop level we get two terms 




ft 



h 2 



l'2 j.fc-5 



2!2!2! (Jfe-5)! 



Lit Ak-4 



3! (k - 4)! 



Finally the last term, i.e the term with the largest number of loop attachments, will simply 
read 



n 



fc-2 



The result looks like 




+ . . . + 



The equation reads 
x = n<f> + Afc ( — — + n 



(k-l)\ 



2! (k - 3)! 



,fc-5 



2 0" fc " 4 



2!2!2! (jfc-5)! 3! (A; -4)! 



+ ... + H 



fc-2 



(F^T)I 



(2) 



For an arbitrary theory a sum over k should be understood. It represents a non-linear 
differential equation for 0, the Schwinger-Dyson (SD) equation, which has been derived by 
the direct application of the Feynman rules. 

In order to be more specific let us consider a theory with only a 3-point and a 4-point 
vertex. Following the abovementioned reasoning a diagrammatic equation for looks like 



x- 




This reads 



x 



4>{x) = - - ^ [ct>{x) 2 + h4>'{x)] - [<P(x) 3 + 3h<j>(x)<j>'(x) + h 2 <f)"(x)} 



(3) 
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This equation generates equations for the Green's functions if the power series of (f>(x) is 
inserted. The first three are 



C\ = — — Ci (A 4 Ci + 3A 3 ) — — C 2 (\aCi + A3) — — A4C3 , 
6/x 2/x 6/i 

1 ft fv^ 

C2 = — — (2A3C1C2 — 2 + A4C1C2) — — (A3C3 + A4C1C3 + A4C2) — — A4C4 , 
2/1 2/j, bfj, 

C3 = — — (2A3C1C3 + 2A3CI + \4C3Cl + 2\^C\C\) 

- ^(3A 4 C 2 C 3 + A3C4 + A4C4C4) - ^A 4 C 5 . (4) 

The SD equation is invariant under certain redefinition of the parameters involved. It is 
not difficult to prove that if \ k , h; x) is a solution, also 

is a solution for any c, a, (5. This scaling property is also a concequence of the fact that 
Afc, h;x)/y/h/fx is a dimensionless function of the scaled variables y = x/^/h/j, and 
gfc = \kh k ^ 2 ~ l I '/i fe ^ 2 - The scaling property can be expressed with the following equations, 
derived, for instance, by differentiating with respect to c and taking c = 1: 

d d , d . d \ , 
ox 011 on oXk J 

r\ r\ r\ \ 

1 + x— + 2/i— + k\ k — 10 = 0. (5) 



dx dfi d\ k 

These equations are equivalent to the usual topological relations that relate the number 
of external lines E, the number of internal lines /, the number of k- vertices V k , and the 
number of loops L, appearing in any diagram, 

kV k = E + 21 V k = I + 1 - L . 

A sum over k should be undesrtood in the general case. 



2.2 Stepping equations 

In the diagrammatic construction, one assumes that every Green's function can be written 
as a sum of diagrams, consisting of vertices connected by lines (propagators). The power 
of in a diagram is equal to the number of propagators, and hence the operation —d/d/i 
on this diagram corresponds to cutting a single propagator in all possible places in that 
diagram. There are two possibilities for the result: the chosen propagator may be either a 
part of a loop, in which case the diagram remains connected when we cut this line, or part 
of the 'tree skeleton', such that cutting it makes the diagram diconnected: 



In the first case the cut diagram remains connected but gains two external lines at the 
price of one loop (i.e. on power of H); in the second place, the cut diagram falls apart into 
two connected diagrams: 

Putting in the correct symmetry factors, we can express this procedure by 

S2: + <Kx)4/(x) + \<j>"{x) = , (6) 

where the second term comes from diagrams that fall apart under cutting, and the third 
one from loops that are cut open. We call Eq.(|6]) the Step-2 equation (S2) since it describes 
a procedure in which the number of external legs is increased in steps of 2. 

Like SD (Eq.(|3|)), the Stepping equation S2 implies relations between various C"s. The 
lowest few of these read 

C*3 — — 7" \C\C2 + ~7\ Ci 



2 / 2 d 

C4 = — — I C1C3 + C 2 + ~qJ^2 

C5 = — — ( CIC4 + 3C2C3 + — — C3 
n \ of! 

C 6 = -^C 1 C 5 + 4C 2 C4 + C 3 2 + |-C^ , (7) 

and so on. Note that S2 is completely independent of the interaction potential, and there- 
fore perforce contains information independent of that contained in the SD. It follows that 
there must be solutions to SD that do not obey S2 and that these solutions cannot be 
represented by Feynman diagrams. 

It is possible to combine SD and S2 in the following manner. Taking the first equation 
in Eq.(||), we express C3 in C2 and Cx, and solve for C^'- 

2A 4 ^ - ~ (6/id + 3A3Q 2 + KCX) 

3A3 + A4C1 

Inserting this into the second equation of Eq.(||), we find a differential equation for G\ 
alone: 

= (?^) 2 - 2h 2 Xl^l (A4C1 + 3A 3 ) 



dfx J d/j, 2 

- n-g± [9A 2 (3A 3 + A4CO + Q/jXld - 36^A3A 4 ] 

- 9/iA 3 C 2 (3A3 + A4C1) + 3h\ 2 4 Cl - 54/i 2 A 3 Ci - 27h\l . (8) 
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By inserting the series expansion 



k>l 



we can then successively determine the coefficients: 
1 



«4 



1 



24/i 5 
1 

"4V 
1 



A 3 , 



(15X1 ~ 16/iA 3 A 4 ) , 

(90A^ - 185/zA^A 4 + 66/i 2 A 3 A^) , 



1152/i 



-(9945A3 7 - 30270^^4 + 24280/i 2 A^A^ - 4352/i 3 A 3 A^) , 



and so on. 

Whereas S2 is independent of the interaction potential, we can also derive stepping 
equations by deleting vertices rather than cutting lines. For example, let us depict all 
possible ways in which a selected f 3 vertex (denoted by a dot) can occur in a connected 
graph: 




Deleting this vertex gives us 



d 



<9A, 



+ 



( ) + 



) + 



or, in terms of (f)(x), the following Step-3 equation (S3): 
S3: 



(9) 



A similar treatment holds for tp A vertices: the possible ways in which such a vertex can 
occur is given by 



and the result of deleting is given by 

d ^ 



<9A 4 




+ 



+ 



+ 



+ 



The corresponding Step-4 equation (S4) is 
S4: + + ^h 2 [200"' + 50V] + \n [0 2 0" + 20(0') 2 ] + = . (10) 
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2.3 The charged scalar field 

Up to now we dealt with diagrammatic construction of zero- dimensional field theories 
involving only one field. As an illustrative extension, we consider a theory with two fields, 
i.e. a complex, or charged, scalar field. The Green's functions are labeled with two integers, 
and the generating function has two expansion parameters x and x. Let us introduce the 
notation 



d-= — 

dx 



d-= — 

dx 



then 



(f)(x,x)= V" ^Cn,m+l • 

n\ ml 

n,m=0 

To write down the SD equation, we introduce two kind of lines, distinguishable by an 
arrow. The generating function is represented by 

(f)(x,x) = ' 4>(x,x) = ~*-(^) ■ 

An incoming external line represents a d, and an outgoing line represents a d. Notice that 

hd(j) = hd(j) = 



We also introduce a four point vertex with two incoming and two outgoing lines, so that 
the SD equation we want to satisfy is given by 



or 



= £__^_« 



2/1 



2/x 



(990 



[11) 
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Notice that incoming and outgoing lines are not equivalent, which is represented by the 
symmetry factors. 

Also for the charged scalar field, we can write down stepping equations. For the first 
one, we use that in the diagrammatic interpretation 

leading to 

A = -hdd(f) - <pd<p - 0d(j) . (12) 

<7/i 

The diagrammatic derivation of the stepping equation involving the derivative with respect 
to A, although equally straightforward, is rather cumbersome, leading to many terms which 
we refrain from listing here. 



9 



3 Solutions to the equations 
3.1 The integral representation 

The SD equation is highly non-linear. Let us consider the general term connected with the 
coupling Afc in Eq.@: 

Qk = 7T~ 7\7 + rtl tl n\ i + ft Jn.n,/, r x, + ^ ^7 , , + • • • + ft fe ~" 



(&-1)! 2!(fc-3)! 2!2!2! (k - 5)! 3! (Jfe — 4)! (k-l)\ ' 
It can obviously be organized such that it can be written as 

ft - S S (1! ).. ai!(2!) ., a2! ...((.- D!).,-^,! w w • • • (^'" 2, )- t - ■ («) 

m=l {a fc _ i; m} 

where stands for the summation with a\, ai, ■ ■ ■ , ak-i running over all positive 

integers under the restrictions that 

ai + 2a 2 + 3a 3 + . . . + {k — l)a&_i = fc — 1 and ai + a 2 + . . . + a^-i = m . 

This sum can be interpreted following the time- honored formula of Fad di Bruno ||: 



f(g(x)) = Y,f im) (9(x)) (n;a 1} ... } a n ){g'(x)r{9"(x)r ■■■{g^(x)Y 



dx n , 

m=0 {a n ;m} 

where 

??.! 



(n;ai,... , a„) 



;i!) ai ai!(2!) ffl 2a 2 !---(n!) a "aj 
The identifications g'(x) = <j)(x) and f (jn \g{x) ) = h~ m f(g(x) ) with the solution 

g(x) = I dx <p(x) , f(g{x)) = R{x) = exp ( - / dx(f)(x)) , (14) 



lead to an equation for R, which, including all possible vertices, reads 

A 



V — — — — rr h k ~ l R( k ~ l) + uhR'ix) -(x- XjRtx) = . (15) 



k=3 

This is a linear equation, and a solution can be represented by an integral 



Rr(x) = J dip exp | — [xip — S(ip)] 



(16) 



where 



1 00 A 



k=3 
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and where T is a contour in the complex 93-plane, such that the difference between the 
values of the integrand in the end-points is zero. This is the well known path integral 
representation, with the action S. 

A remark is in order. Although the SD equations resulted from a a purely diagrammatic 
construction, their solutions, expressed through the path integral representation, include 
non-perturbative ones that cannot be realized in a weak coupling expansion, as we will see 
below. 

Secondly, we note that Eq.flToT) can be used to write the original SD equation for (f> 
compactly as 

i / o \ k-2 



k>3 



For a general interacting theory, differantiating R-p with respect to A& in Eq.(|TB|), the 
stepping equation in terms of <fi can be rewritten as 

d(j) 1 d f d N A ' 



d\k k\ dx \ dx 

and in case only a fc-vertex and the tadpole Ai is present, combining with the SD a simpler 
form is obtained 

lr = -ir ( - Ai) 0' + - 2^00' - w') • (19) 

Moreover for a charged scalar field the stepping equation in terms of <ft can be written in 
a compact form 

^ = - l d{4> + hdf{ ( t) + hd)(t) . (20) 

OA 4 

Finally, the linear SD equation for (f 3 + <^ 4 -theory becomes simply 

l\ 4 h 3 R"'(x) + jA 3 ^ 2 i? // (x) + iihJUix) - xi?(x) = , (21) 
6 2 

and we see that R(x) admits 3 linearly independent solutions (2 if A4 = 0). Hence <f>(x) 
has a 2-parameter family of solutions (a 1-parameter family if A4 = 0). In the sequel we 
will show how to get exact explicit solutions for a number of scalar theories. 

3.2 Results for pure <£> 3 -theory 

In this section we derive results for the pure v? 3 -theory, with action 

S{<p) = \ W 2 + ^ 3 . (22) 
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This theory is interesting because as we will see the solution for the generating function 
can be expressed directly in terms of known special functions. Defining 

x xVh 
y = ~n= > ? 



/ K/2 ' 6/i 3 / 2 
the SD equation becomes 

3^R"(y) + R'(y)-yR(y) = (23) 
which admits the following general solution 

R(y) = e-^[ Cl M(t) + c 2 Bi(t)} (24) 

where 

t = m- 1/3 i^ + y 



12£ 

Ai and Bi are the Airy functions (cf. [[J). The solution for the generating function of 
connected Green's functions is given by 

V A* V 6£ Ai(t) + ABi(t) y 1 ; 

with t = t(y = 0). The constant K is not determined by the SD equation: in fact it could 
have been even a function of £. 

For solutions that admit a diagrammatic representation extra information can be ob- 
tained by combining SD and stepping equations. For instance, the scaling and stepping 
equations of the previous section result to a A that is independent of £. Moreover by 
combining SD and S2 an equation involving only C\. 

2/j 2 Ci + AyuC? + hX 2 -^-C 1 + h\ = , (26) 

can be obtained. The series of substitutions 

v = ^ , Ci = ~f{v) /(„) = , fe(«;) = ^/V^M , 

/i d 2/i z 6v vk[w) 

leads to the Bessel equation 

w 2 ip"(w) + wip'(w) - (w 2 + ^ = • 

The special solution choice ip(w) = K 1 / 3 (w) gives the following tadpole and its asymptotic 
expansion: 

2 fK 2/3 (w) 



/(*) = " 1F^K + 1 
v \Ki/ 3 (w) 



a h\ 2 ( 5 15 2 1105 3 1695 A 



C 2^ - — - 1 - + -V - + + ■ ■ • . (27) 

A 2fi 2 V 4 4 64 16 ' 
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This tadpole, therefore, has a non-perturbative contribution. The more generic choice 
ip = kil 1/3 (w) + k 2 I- 1 / 3 (w), with k x ^ -k 2 , gives 



2 fkj^/^w) + k 2 I 2/3 (w) 
f(v) = — 



v \k\h/z{w) + k 2 I_i /3 (w) 
h\ 2 / 5 15 2 1105 3 1695 4 \ 

which is the standard perturbative result j|. The coefficients fci and fc 2 drop out for the 
perturbative expansion: they simply account for non-perturbative contributions that are 
not computable perturbatively! 

A remark is in order here: although all the terms in the perturbative series for C% have 
strictly the same complex phase and according to the traditional wisdom the series is not 
Borel summable, the exact result is well defined, indicating that a suitable generalization 
of the Borel transform will produce the right answer ||. 

Another interesting aspect is the large n behavior of the Green's functions, where n 
refers to the number of external legs, a problem that is traditionally seen as relevant to 
the unitarity of the S matrix ||. This can be traced from the analytical structure of the 
solution for the generating functions in the complex x-plane. As is evident from the fact 
that the solution, Eq.(plD, for the generating function of all connected and disconnected 
graphs is an entire function, the corresponding Green's function Z n grows slower than n!; 
in fact it grows like (n!) 2 / 3 . On the other hand the C n , the connected graphs, exhibit 
a factorial growth, since their generating function <j>(x) posesses poles at finite complex 
values of x. 

3.3 Results for pure <£> 4 -theory 

In this section we derive the lowest Green's functions for the pure f A theory, with action 

S(<p) = \w 2 + ^Ay, 4 . (29) 



Defining 



x 



a T/hjl * 24/i 2 

we get for the SD equation 

^R"\y) + R'(y)-yR(y) = . (30) 
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There are three solutions, which can be represented as follows: 

oo 2n 

Ri(y) = E V -y (320~ n/2 u(n ; (so- 1/2 ) 

J 2 " /90 ^ n / 2 V("i (8Q- 1/2 ) 



n! 

n=0 

oo 



n=0 ^ 2 ' 

oo 2n+l 

(2n + l)! 



n=0 

where U(f; x) and V(v; a;) are the parabolic cylinder functions, and H n is the n th Hermite 
polynomial (cf. 0). The general solution is a linear combination with arbitrary coefficients. 
As we can immediately see contrary to what is argued in many standard textbooks the 
odd Green's functions do not necessarily vanish. 

On the other hand on can study the S2 equation as well. For this we have to distinguish 
two possible cases: the 'standard' one, with C\ and the higher odd Green's functions 
vanishing, and the case where C\ ^ 0. 

Let us consider the first case with a zero tadpole. In this case we cannot, of course, 
directly use the results derived above, since these deal with C\. The S2 becomes somewhat 
simpler, and in particular 

Ci = ~\ (5 + 
On dimensional grounds we see that we can write 

C 2 = ~(3(v) , v = \ , 

where v is dimensionless. Inserting all this into the first nonzero term (that with x 1 ) in 
SD, we find the following equation for (5: 

Av 2 (3\v) + v(3(v) 2 + (2v + 6)p(v) -6 = 0. 

The substitutions 

i3(v) = Av^- , g(v) = v-^e^iw) , w = f , (32) 
g{v) Av 

lead then to 

ii) 2 f(ff) + wi//(w) - (w 2 + ij) (w) = , (33) 

which has the modified Bessel functions for its solutions. The general solution can always 
be written as 
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It is instructive to consider the perturbative form of these results, that is, the limit where 
h becomes infinitesimally small, or w goes to infinity. Since Iiu and I-1/4 have the same 
asymptotic expansion, a generic choice of kip will lead to a single perturbative expan- 
sion. The single exception is the choice h\ = —k 2 which leads to ip(w) oc Ki/±(w), with 
asymptotic expansion 

K V*W ,\ 1 1 2 2 11 3 34 4 



' y ' v \K 1/4 (w) J 2 3 8 9 

This is the standard perturbative expansion, in which the propagator starts with ^, and 
has loop corrections in powers of H: 

The alternating signs are of course due to the fact that the Feynman rules prescribe a 
factor —A for each vertex in our Euclidean model. The asymptotic expansion in all other 
cases is equal to that for the choice k 2 — 0, for which we find 

3 //-3/4W A 6 1 1 2 2 11 3 34 4 

P(v) = — T 9 \ + 1 + l + -v + -v 2 + -v 3 + -v 4 + --- , 

v \ Ii/i\w) J v 2 3 8 9 

which gives a nonstandard expansion: 

Xh n\ 2/i 2 3/i 



C = d? = -^ + i|l + a + ^f- + -) ■ (35) 



Note the occurrence of a 'non-perturbative' term 1/A here: the rest of the expansion has an 
apparent opposite sign of the coupling constant. An other way to look at this solution is by 
examining the saddle point equation, 5S/ 5<p = x: the abovementioned solution corresponds 
to the saddle point (j) c = a/— 6^/A + 0(x). 

In the case C\ 7^ we can write, again on dimensional grounds, 

Ci = a(v)^ , C 2 = , 

with v as before. The first term (with x°) in SD now gives us a relation between a and f3: 

P(v) = — ^ ((6 - v)a(v) + 4v 2 a'(v)) , 
va[v) 

and then the second term (x 1 ) gives 

16v 2 a(v)a"(v) - 32v 2 a'(v) 2 + (32v - 2A)a(v)a'(v) - 3a(v) 2 = . 
Using w as before, we may now substitute 

ip(w) 
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to find that ip(w) again obeys the Bessel equation, Eq. (|33|) . For the asymptotic expansions, 
again two distinct choices are possible. First, the choice 

1 

ip(w) = -K 1/4 (w) 



gives 



pe"yfr 3 (K z/A {w) ^ p 2 e 2w 



K 1/4 (w) ' v \K 1/4 (w) J K 1/4 (w 

and the following asymptotic forms for C 4 ^'- 



V 2nX fin 



The alternative choice, for which we may take 

i>{w) = -h /4 (w) 
p 



leads to 



pe w J6 3 ( I-y 4 (w) \ p 2 e 2w 

a ( v ) = 7 — 7 \ > = — - 7 , , + 1 - 



Ji/ 4 (w) u V h/i(w) J h/ 4 (w 



and 



/3tt/7 (2) 2vrp 2 w / 1 



IV _ Z 3 ^ _ . . . ] (37) 
96 384 ' 



In contrast to the zero-tadpole case, there remains an arbitrary parameter in these solu- 
tions, p: it reflects the presence of the 'non-perturbative' tadpole-like contribution and has 
to be determined by additional requirements. 

3.4 Results for cp 3 + cp 4 - theory 

For the general zero-dimensional (/2 3 + y9 4 -theory, the action is given by 

%>) = \w 2 + + y a x ^ 4 ■ ( 38 ) 

In the dimensionless variables 

x A3 h \ 4 h 

v = —f^ > 93 = —\\-i 9a = — , 



the SD equation becomes 



\g 4 R"\y) + \g,R'\y) + R'(y) - yR(y) = . (39) 
6 2 
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To solve this equation, let 

R( y ) = e- y93 ' 9A F{y) 

Then F satisfies the equation 



where 



\gtF"'{y) + aF'(y) - (y + (3)F{y) = , (40) 
o 



2^4 9i \ 3^4 



Finally, changing variables 
Eq.pDD becomes 

4£F"'(r)) + F'(r}) - rjF(r}) = . (41) 



Eq.(f4"l~|) is exactly Eq.(^) of the pure </9 4 -theory, so that the solutions here are those of 
([31]) with £ as given above and y replaced by rj. 

3.5 Results for the charged scalar field 

For the complex scalar field, the path integral solution is given by 



R(x, x) — J dipdip exp | — [xcp + xip — S(<p, (f))j 



\ _ 2 

S(<p, <p) = (i<p<p + - (<p<p) . (42) 



Due to charge conservation (0(2)-symmetry) one can easily show that R only depends on 
the modulus Ob Ob . SO that it satisfies the following equation 

CR'"(C) + 2R"(C) + a 2 {R'(C)-R(C)} = o , C=~, a = /i (^) V2 • (43) 

The third order equation can be solved by power series expansion in ( and two of its 
solutions are given by 

*<o=E^gf)"B.(^) • «.(0 = E^«(» + 1,-), <«) 

n=0 \ v / \ v / n=0 

where H n stands for the n th order Hermite polynomial and U(f, x) is the parabolic cylinder 
function. The third solution can be found by standard procedures but the actual result is 
rather cumbersome and we refrain from giving it explicitly. 
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Where R\ is a purely non-perturbative solution, i? 2 has an asymptotic series expansion 
which leads to the normal perturbation series for the connected Green's functions: for 
instance the two point function is given by 

C, 1= l --^ + \^ + ( m - (45) 

Although the standard integral representation of R is given by the differential equation 
in the variable ( leads to another peculiar single contour integral representation 

*0 = jU«p(»-ln*-J + 5j) , 
where the contour F is from infinity to infinity such that the integral is convergent. 



3.6 Contours in the integral representation 

The integral representation of the solutions for the pure v? 3 -theory, with ip = a/ /J,/h(p, can 
be written as 



R(y,0 = K jf#exp 



^_^3 +# ; (46) 



where K is a constant which can depend on £. In case the moduli of the endpoints of the 
contour V are taken to infinity the standard path-integral representation is recovered. In 
fact in this case the substitution 

u = (30~ 1/3 ^ - ^ 

leads to 

R(v,0 = Kexp( — | / du exr> ( — u 3 + tu] . 

This integral can now be expressed in terms of the Airy functions 

J du exp (^-^ u3 + tuj = 2nuj j Ai(ttu j ) , 

where u = e i27T/3 , j = 0, 1, 2 and the contours Tj are depicted in Fig.|l|. Note that 

^Jki(tJ) =0 . 

3 

For a pure y9 4 -theory, similar considerations allow us to express the functions Rj defined 
in Eq.(0), j = 1, 2, 3 as follows 

Ri = i(20 1/4 exp ( ] / di) exp ( ~ip 2 - ^ + yip 



7T V J-oc V 2 



R 2 = -A(20 1/4 exp f-^M J™ # exp f-~V> 2 - & 4 + yij\ (47) 
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Figure 1: Contours in the complex w-plane for the Airy functions. 



whereas for R 3 

{Jo " + I tQO ) # ^ ("^ ~ ^ 4 + # ) 

= ^(20- 1/4 exp (ifc -M 2 ) + |exp (J*)" 1 ' 4 ^ . (48) 

Let us have a closer look at the various possible contours in the case of general tp 3 +Lp 4 - 
theory. Let us denote the various objects in the action as complex numbers: 

up = \ v \ e luJ , A 3 = I A 3 |e^ 3 , A 4 = |A 4 |e ir ? 4 . 

For simplicity and without loss of generality, we may keep \x real and positive. The direction 
in the <p plane where the term A 4 y? 4 goes to positive infinity as \ip\ ^00 are given by 

ueQAk) , nAk) = (k---- 7 ^ , k- + --^) , k = 0,1, 2, 3 . 

4V ) > 4 V ; ^ 2 g 4 ' 2 8 4 / 
Similarly 'allowed' directions for the A39? 3 term are 

Finally, the /i<£> 2 term goes to positive infinity for 

e fi 2 (/c) , n 2 (A;) = (kir - j , A;7T + ^ ,fc = 0,l. 

By inspection of these endpoints, already statements can be made about the (non)pertur- 
bative character of the theory corresponding to a given contour. To illustrate this, let us 
consider a pure y? 4 -theory, i.e. with A 3 = 0. Let the contour start at some ipx, chosen 
at infinity with argument u>i, and end at some ip 2 , also at infinity in some direction with 
argument u>2- These values each have to be in some interval f2 4 : let uj\ be in f2 4 (ni), and uj 2 
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Figure 2: The regions in the complex </>plane which correspond with Q^k), Qs(k) and 
f2 2 (/c) with rjj = —joij. 



in fi(n 2 ). We can sufficiently specify the contour by giving n\ and n 2 so that for instance 
the contour r 20 for 774 = denotes the standard <^ 4 -theory, where we may take the real line 
for T, start at ip — — oo and end at tp = +00 (interchange of the endpoints corresponds 
to replacing R by —R and hence does not influence (f>(x)). In total, there are six contours 
that give a viable (//-theory: r i, T 12 , T 23 , r 30 , T 02 and T 13 . Note that these are related 
to each other by phase shifts: in fact, 



r 3 o = ri 2 (r? 4 -> 774 + 27r) , r 23 = 1^(774 -> 774 + 47r) , 
r i2 = ri 2 (r] 4 -> 774 + 6tt) , r 13 = r 02 (r] 4 -> ^4 + 2tt) . 

Therefore, only T 02 and r i, say, give really different theories, all other cases being obtain- 
able by an appropriate shift in 774. All contours, as stated, corresponds to viable theories 
as long as A4 is non- vanishing, but when we let | A4 1 — > there are two possibilities. It 
may happen that ^4(711) overlaps with one of the fl 2 segments, and ^4(n 2 ) with the other 
Q 2 segment. In that case, the limiting theory is equal to the free theory, and the limit 
I A4 1 — >• is smooth: we may call this the perturbative limit. In the other case the limit 
is not smooth, and the path integral R will diverge as |A 4 | — > 0: we call this the non- 
perturbative limit. Clearly, the limiting behavior depends on the argument 774: for the 
contour T 02 (the 'standard one') one has 



--rr< m <-rr 

3 5 
2 n<m< 2 n ' 

and for the other contour Tqi: 



-7T < 77 4 < -TT 
2 ' 2 



5 11 

2 n<r]i< Y n 

11 13 

Y n< V4 < Y n 



9 11 

2 ' 2 



11 
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-7T < 7]^ < -TT , —11 < ?7 4 < yvr : 



perturbative , 
non-perturbative , 

perturbative , 
non-perturbative , 
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Figure 3: The shaded areas correspond to combinations of 774 and 773 (in units of 7r) for 
which the limit A 4 — > (with A 3 fixed) is perturbative, for contour r i (left plot) and 
contour T 2 (right). 

For the pure y2 3 -theory, there are of course three contours, related to each other: r 20 = 
roi(^3 — ► V3 + 27r), Ti2 = r i(?73 — > 773 + 471"). For the limiting theory we find, for contour 
I01: 

5 1 7 13 

— -7r < 773 < -7r , -7T < ^3 < — vr : perturbative , 

1 7 13 19 

-7T < r/3 < -7T , — 7T < r/3 < — 7r : non-perturbative . 

In a theory with both y) 3 and (f A couplings, things become more interesting. Of course, 
as long as A 4 is nonzero, we are allowed to let A 3 go to zero without jeopardizing the 
perturbativity. On the other hand, we can only let A4 vanish with fixed A3 if the selected 
f^4 and Q3 intervals overlap. We give in Fig.|3] and Fig.^ the values of 774 and r/3 that 
correspond to a perturbative A4 — ► limit, Finally, we may study the combined limits 
A 4 — > followed by A3 — > 0. The regions of perturbativity are given in Fig||. Clearly, these 
are more restricted since we are in this case requiring a common overlap of Q2, ^3 and Q4. 
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Figure 4: The shaded areas correspond to combinations of 774 and i] 3 (in units of ir) for 
which the limit A 4 — > followed by A 3 — > is perturbative, for contour r 01 (left) and r 02 . 
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4 Renormalization 



Even though zero- dimensional field theories have no infinities, we may still consider the 
effects of renormalization, which here take a graph-theoretical significance. Renormalizing 
the field so that the exact propagator is l//i, and the coupling constant so that the proper 
vertices assume their tree-order form, we are counting Green's functions without self-energy 
and vertex insertions, that is, we are counting the skeleton diagrams of the theory. We 
shall restrict ourselves to theories that are known to be perturbatively renormalizable in 
the usual four-dimensional case. 

4.1 Renormalization of pure <£> 3 -theory 

In this case renormalization proceeds as usual with the introduction of the tadpole (z\), the 
mass (Z2) and the vertex (2:3) counter terms, as well as the corresponding renormalization 
conditions that imply the dependence of these counter terms on the renormalized coupling 
constant. The renormalized action can be written as (/i = 1, % = 1) 

1 1 

S = -z 2 <£ 2 + -^gz 3 ^ 3 + z X ip , (49) 

2 



and the SD equation takes the form 

z 2 <p=(x-z l )-^(<f ) 2 + <p') , G 3 = gz 3 . (50) 
Moreover using Eq.(|T9"[) we get, 

3G 3^ = W + 2^200' - - (z - Z X )<jf , (51) 



whereas Eq.(pl) becomes, 



# = -00' - V 

dz 2 vv r 



(52) 



The renormalization conditions that have to be applied are 

Condition 4.1.1. No tadpoles, i.e. 0(x = 0) = ; 
Condition 4.1.2. propagator = 0'(O) = 1 ; 
Condition 4.1.3. vertex = 0"(O) = —g. 

Application of these conditions to the SD equation and its derivative leads to the equations 

z\ = » z 2 = l + ^g 2 z 3 . (53) 
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So if we know z 3 as function of g, we know z\ and z 2 as function of g, and we can consider 
to be a function of g and x only. Its derivative w.r.t. g can, using d/dz\ = —d/dx, be 
written as 

^ = -^ 1 + ^ 2 + ^ Gs (54) 

where a dot denotes differentiation w.r.t. g. Because is a function of x and g only, the 
l.h.s. is zero in x = by Condition [4.1.1| , and evaluation of the r.h.s. leads to the equation 



gdz 3 2z 3 - 2(1 + g 1 )z, 



2^2 
3 



2dg -4 + (4 + g*)z 3 
It is straightforward to derive that for g = the perturbative counter terms read 

2i(0) = 0, z 2 (0) = l, 2a(0) = 1. 



(55) 



Eq.(^5]) is an Abel equation of the second kind J7J. The perturbative solution, satisfying 
the above initial condition, is 

zs(g) = l-9 2 - 1 -9 A - ^ ~ 2V - 5 -fg w + ■ • • , (56) 

an expansion previously given by Cvitanovic et al. ||. 

Having, however, solved the SD equation for ^-theory (Section |3.2| ), we also have an 
exact, albeit implicit, solution of Eq.(^), making use of Condition \l.l.2\ 



1/3 

\gz 3 ) gz 3 



[ Cl Ai'(t ) + c 2 Bi'(t )] (—) - — [ciAi(to) + c 2 Bi(t )] = , (57) 



where 



/ 2 \ 1/3 A z\ \ z\ ( g 2 z\ 

to -{-gV 3 ) \2 9Z * + 2^ 3 )-WWA l + ^¥ 

and Ai and Bi are the two independent solutions of the Airy equation f"(t) = tf(t). The 
meaning of this equation Eq.([57]) is that for a given g and by using Eq.(|53l) as well as the 
functional form of Ai and Bi we can determine z 3 . To show that Eq. ([57D is an implicit 
solution of Eq.(|55|), let 

rn/ , (2g 2 zl) 1 / 3 

implying toF(g) 2 — 1 = (gz 3 /z2) 2 , and differentiate Eq . (|57|) with respect to g to get 
(F'(g) - ^ [ Cl Ai'(t ) + c 2 Bi'(t )] + F(g)[ Cl Ai"(t ) + c 2 Bi"(t )]^ = . 
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2 4 6 8 10 -0.2 0.2 0.4 0.6 0.8 1 1.2 

Figure 5: z% as function of Re g with Img = for pure y? 3 -theory. 




2 4 6 8 10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 6: g as function of Re gz 3 with lmgz3 = for pure </> 3 -theory. 



Using the Airy equation and Eq . (|57|) , we get 



Explicitly, this says 



F'{g) + ^(toHg) 2 - 1) = 



2z 2 ^ - 3gz 3 ^ - 2z 3 ^ = 
dg dg dg 



By using ([53]) , one easily sees that the above equation is an equivalent form of Eq. ( 

In Fig.[5] we present the results of a numerical calculation of Z3 for the Tio-contour as 
function of g, as described in the Appendix. The left graph shows the real and imaginary 
part of z 3 as function of real and positive values of g. Notice that z 3 (0) = 1 as demanded, 
and that the imaginary part does not stay zero for real g. This is, of course, an artifact 
of the definition of the path integral over a complex contour, which is the Fio-contour for 
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y9 3 -theory in this case (Fig.|2|). The right graph combines the real and imaginary part in 
one curve in the complex z 3 -plane. 

Fig.^ shows what happens if we let g run with real and positive values of gz%, so that 
the actual coupling constant is real and positive. 

4.2 Renormalization of pure <£> 4 -theory 

In the case of <^ 4 -theory, the renormalized action is given by (// = % = 1) 

S = \ z 2^ 2 + J^z^ ■ (58) 

The SD equation becomes 

^ 2 = x-^(0 3 + 300' + 0") , G 4 = gz 4 , (59) 
o 

and the stepping equation Eq . (|T^) , leads to 

4G 4 = z 2 <\>" + 2z 2 #' -<f>-x<f>' , (60) 

whereas Eq.(^) assumes the form of Eq.([52|). The renormalization conditions require that 

Condition 4.2.1. <p(x = 0) = 0; 
Condition 4.2.2. 0'(O) = 1; 
Condition 4.2.3. 0"'(O) = -g, 

and application to the SD equation leads to the relation 

z 2 = 1- g)gz 4 . (61) 

6 

As in the case of </2 3 -theory, can be considered to be a function of x and g only, and its 
derivative w.r.t. g can be written as 

d(j) _d<P d<j) ■ 

dg~-dz- 2 Z2 + dcf* ■ (62) 



The l.h.s. is zero in x = by Condition |4.2.1| , and evaluation of the r.h.s. leads to 

dz A -6z 4 + (6 - 9g + 3g 2 )zj 
dg 6# - g(6 -5g + g 2 )z 4 

another Abel equation of the second kind. The perturbative solution is given by 



(63) 



, \ 1 3 3 2 11 3 45 4 499 5 
zM = l + - 2 9+-/ + Y c f--t ? + f + ... . (64) 
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We want to remark at this point that the statement by Cvitanovic et al. that, in 
the case of <^ 3 -theory, the coefficients of the series expansion of g — gz 3 count connected 
three-point diagrams with no self-energy or vertex insertions cannot be carried foreward 
to <£> 4 -theory: the coefficients of the series expansion of g — gz^ do not count connected 
four-point diagrams with no self-energy or (four-point) vertex insertions. There are, for 
example, no such diagrams with three vertices. 

To find the exact implicit solution of Eq . (|63|) , we apply Condition |4.2.2| to the solution 



(|3T|) of (f -theory, resulting in 



CitU(l;t) c 2 tV(l;t) TT ,_ c 2 V(0;t) f 3z 



:! - Cl U(0; t) + > * = ( ) • ( 65 ) 



2\ 1/2 



Z2 z 2 T(l) x v w r(±) \gz 4 

Letting 

Fi(t) = Cl U(l;t) - ^1 and F (t) = Cl U(0;t) + ^1 
the above equation becomes 

tF 1 (t) = z 2 F (t) . 

Using the properties of the parabolic functions we can easily show that 

Fl(t) = (| - l) F (t) , F '(t) = ~ (t + I) F (t) , 

so that differentiation of Eq. ([4.2.3| ) leads to 

dt z 2 „ / s dt ( z 2 \ . . dz 2 i \ z 2 dt / z 2 \ „ , . 

*f F ° w + % a - o F » (i > = -i F » (i> - f * (« + f) F ° (i » ■ 



This equation can be written as 

d t dt (\ t 2 g-3 \ _ 
dgz 2 dg \2 z\ 6 4 / 

where we used relation fl61|). Finally, since t/z 2 = a/3/ ^4 by definition of t, it is easily 
seen that the above equation becomes Eq.flBBp. 

In Fig. [I] we show the results of the numerical calculation of z^(g), as described in the 
Appendix. We used the r 2 o-contour for ^-theory (Fig.^j) . Starting at 24(0) = 1, z±{g) 
stays real and positive for real and positive values of g, as expected. Moreover z 2 exhibits 
a zero, whose position g+ can be calculated analytically and is given by 

^ = 3-M^||| -0.81155 
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Figure 7: z± as function of g (left) and z 2 as function of g (right) 




Figure 8: C 6 as function of g (left) and logz 4 as function of — log(2 — g) (right). 
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-200 



Figure 9: The complex g-plane and the complex 2 4 -plane. g goes around twice, anti- 
clockwise and starting on the real axis on the left side of g = 2. 



At this point the theory becomes 'massless', in the sense that the bare mass becomes zero, 
yet the Green's functions do not exhibit singular behavior. In fact let us consider the 
6-point function as an example. It can be explicitly calculated and it reads 

C 6 = 6z^ -6 + 9g + g 2 . 

It is easy to see that the expansion around g = reproduces the known perturbative series. 
Moreover, the left graph Fig.|8| presents Cq as a function of g. 

We also see that 24 increases with increasing g, and explodes if g approaches 2. The 
graph on the right of Fig.|] suggests that the there is a simple pole at g = 2. In fact, 
substitution of a Laurent series around g = 2 in Eq. ( |63| ) results in a solution with a simple 
pole: 



24(g) 



6 



+ 6-12(0-2) + 54(s- 2) 



399( 



2) 3 + 3948(#-2) 4 - 



(67) 



One can ask the question whether this series expansion corresponds to a solution with 
2; 4 (0) = 1, that is, the perturbative solution. In order to get the perturbative solution from 
the implicit solution fl65|) , in combination with Eq. (|6T|) , we should take the constants C\ 
and C2 such that the limit of t — > 00 exists. Using the properties of the parabolic functions 
and their asymptotic expansions, we find that the perturbative solution has to satisfy 



z 2 =t 2 



l+2\ 



2 sin utt 



[I- 



1+2 



if Ret < 



hat 2 )} ifRet>0, 



(6£ 



together with Eq . (plf) . For g close to, but smaller than, g = 2 we see that z 2 < 0, so that 
Ret < 0, and it is easy to see that the solution in this case has a simple pole at g = 2. 
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However, the coefficients for large powers in the series expansion seem to behave as (2n+2)!, 
so that the series has radius of convergence equal to zero, and the numerical solution of 
a curve around g = 2 in the complex g-plane reveals that there is a branch point (Fig|J). 
In any case, for g —>■ 2~ the bare coupling becomes strong and the bare mass squared 
large and negative whereas the connected Green's functions are still perfectly calculable; 
for instance Ce(g = 2) = 16. 

4.3 Renormalization of ip 3 + <£> 4 -theory 

The renormalization of the 9? 3 -|-<^ 4 -theory is more involved, but straightforward. The action 
is given by 



S = \ Z W 2 + ^[ G 3^ 3 + \\ G ^ + Z W 

and the SD equation assumes the form 



G3 — g-iZ-i , G4 — #42:4 



z 2 <p = (X - Z\. 



+ 3#' + <P") 



The stepping equations read 



d(j) 

da, 



dG 4 



dz 2 



vr^ 4 



6 r 
1 



24^ 6 



and the renormalization conditions are now 

Condition 4.3.1. <f>{x = 0) = 0; 

Condition 4.3.2. 0'(O) = 1; 

Condition 4.3.3. 0"(O) = -g 3 ; 

Condition 4.3.4. 0"'(O) = 3#J - # 4 - 



(69) 



Combining these conditions with the SD equation one easily gets 



—q^Ga — — G3 
2 2 



^2 = 1- -(3flr| - 5-4 + 3)G 4 + -5-3^3 
6 2 



(70) 



so that <j) becomes a function of g$, g 4 and x only, leading to the the four equations: 



dgi 



x=0 



+ 



x=0 d 9i dG 



r=0 9 9i 



d<j> 



dG, 



dG A 

x=0 9 9i 
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and 



x=0 



dz 2 d<fi' 
__ dgi dG 3 



dGo 



+ 



__ dgi dG 4 



dG 4 
=o d 9i 







with % — 3, 4. The coefficients at x = can be inferred form the stepping equations. This 
way we have a system of four equations involving the partial derivatives of the functions 
G 3 (g 3 , g 4 ) and G 4 (g 3 , g 4 ) with respect to g 3 and g 4 . Notice that the equations are linear with 
respect to the four partial derivatives but higly non-linear with respect to the functions 
G 3 (g 3 ,g 4 ) and G 4 (g 3 , g 4 ). They can be solved perturbatively with the result 



3 3 3 

= 93 ~ 9l h + g 9s 94 h + 4 gt ft 2 - 6 g 3 g 4 h 2 + - g 3 g\ h 2 - 4 g\ ft 3 - - g\ g 4 ft 3 

+ ^ //;; 9l I' '' • y//3 <?1 * 3 • T^lr-^- g\ g 4 h* + 81 ^ ^ /, 1 - ^ ^ <?t * 4 



• 927 ^/i 6 =^g-l,flir 



+ 



G A 



45 4 + 4 11*5 807 9 

-g 3 g\fti i + 47g 1 3 1 h 5 - — 
1785 3 4 ^ 5 499 5 ^ 5 

g 4 + 3glh-6glg 4 h + ^glh-6glh 2 + 5g%g 4 h 2 + ^g 2 gjh 2 + ^g 3 h 2 



+ 9g 8 3 h 3 -A3g b 3 g 4 h 



6 



151 



gt g\ h 3 - 39 g 2 gt ft 6 + ^ <? 4 4 ft a + 33 ^ 3 iU ft 



,2 „3 *3 



11 



4*3 



,10 *4 



I . oo , fi .2 ?- I 1485 „1 3 t I 



324^^^ + 834^^^- 
1029 



1585 2 4 4 45 5 4 



2 

L2f.o if.m .10,, /■! , 27525 8 9*5 i 7 non „6 „3 *5 , 68595 , - 



+ —9S 1 * -4610^ 1U 5 4 ti>+ , 

10705 2 5i5 499 6i5 
- —g— 03 04 ft 5 + -33" 04^ 



^^- 17020^^ + 



4 4 f-0 

#3 04^ 



where the ft dependence has been restored for convenience. 

In the limit g 3 — > 0, also G 3 — > 0, and the equations reduce to 



dG 4 

dg 3 
dG 3 

dg 3 





CU 
#4 



9G4 

<%4 

dg 4 



2(2 -g 4 )G 



%4 - (6 - 5^4 + gf)G 4 
. 



Note that (Lt 4 (0, g 4 ) can be identified as g 4 z 4 (g 4 ) where z 4 is the same function is as in pure 
y? 4 -theory. Another interesting result is that the term linear in g 3 in the expansion of G 3 
is given by 



G 3 {g 3 ,g 4 ) = 53^4(54) + 0(g 2 3 ) 
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4.4 Renormalization of the charged scalar field 

In the case of the charged scalar field we consider the integral representation 



R(x, x) = \ ^ I d(pd(p exp < — : 



h J r T { h 

which satisfies the SD equation 



/J,Z2(p(p + ^-^4 (<£><£>) 2 — Xif — X(p 



implying 



dR ld(gz 4 ) d 4 R dz 2 d 2 R 



Condition 4.4.2. ®} A x = x = 0) = -^r , implying R"(Q) = ( 1 - -) R(Q) 
oxox 2 [T V 2/ 



(71) 



CR"\Q + 2R"{Q + — [z 2 R\O-R(()] = , C = ^f, 9=™- (72) 
In the dimensionless variables 

Eq.([7TD becomes 

R(u, u) = I dipdip exp (—Z2ipip — -gz^ipip)' 2 + nip + uip 



dg 4 dg du 2 dv? dg dudu 
or, in terms of the (^-variable 

I? = _ I fcl (fR»» + 4(R"' + 2R") - ^((R" + R') • (73) 

Now, the generating function of the connected Green's functions is given by 

, / \ *. 9 . _ . A . x R'(C) 

= in jkch--^ , 

and the renormalization conditions are 

Condition 4.4.1. ^-(x = x = 0) = - , implying R'(0) = R(0) ; 



By combining equations Eq.([72"D and Eq.(^) and the renormalization conditions we get 



- 2 = I ( 1 - I J gz 4 



32 



and 



dz 4 2z 4 - (3g 2 -5g + 2)z] 



' (74) 



dg -2g + g(g 2 -3g + 2)z 4 ' 

with perturbative expansion 

. , 1 5 9 2 49 3 271 4 5025 5 

zAq) = 1 + -q + -q H o 3 g H g H . (75) 

4 V y; 2 4 8 16 32 y ' 



To get an exact implicit solution of Eq . ([74]) , we go back to Eq . (|72|) and change variables 
to get 



r]R"'(v) + 2R"(V) + a 2 [R'(v) ~ R(V)} = , V = — , « = W — z 2 . (76) 

2 2 V 5^4 



This equation has exactly the form of (|43|) , and the perturbative solution is given by 

^) = E^rU(n + i;a) ■ 

n=0 



Condition [4.4.1| implies the implicit exact solution of the form 

aU(|;a) = z 2 U(±;a) , (77) 

where, of course, z 2 = 1 — (1 — g/2)gz 4 . 

To show that Eq.flT7|) is indeed an implicit solution of Eq. (|T4"[) , differentiate Eq. ([T7|) 
with respect to g: 



J; a) + aU'(§; a)]^ = ^U(|; a) + z 2 U% a)^ 



and using parabolic cylinder functions properties together with Eq. (|77|) to get 

dz 2 (z 2 z\ \ da 

— = Yaz 2 ^ a h- , (78) 

dg \ a a J dg 

with a = ^J^2 z 2 . It is a straight forward calculation to show that this is indeed Eq . (|74]) . 

In the following we present a derivation of the initial condition for z 4 (g = 0). Using the 
path integral expression of Eq.([H]), we find the SD equation 

liz 2 hdR + ^-h 3 dd 2 R-xR = . 

The generating function = hd In R of the connected diagrams satisfies 

0(0, 0) = , (00) (0, 0) = , 0(0, 0) = , (<90)(O, 0) = . 
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as well as the SD equation 



''' A : 1 ' 7 12 + 2ft0«90 + 0<90 + h 2 3d(j)) 



fiz 2 2/iz 2 

with the renormalization conditions |4.4.1| and |4.4.2| . These should hold for any value of ft, 
and for ft = 0, the SD equation becomes 



x 



Xz 4 {0) - 



fj,z 2 {0) 2/i2: 2 (0) 
from which we derive for the perturbative solution that 

1 



(<90o)(O,O) 
(<9d 2 0o)(O,O) 



/j,z 2 (0) 
Az 4 (0) 



I? 



z 2 (0) = 1 

Z 4 (0) = 1 



Notice that the value ft = is directly related to g = since (7 is proportional to ft. 
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5 Summary 



In this paper we studied several aspects of zero- dimensional field theories. In the first 
place we derived a set of diagrammatic equations, including the well known Schwinger- 
Dyson equations as well as a set of 'stepping' equations generalizing some previous results. 
Then we showed how to solve these equations exactly in terms of known functions and we 
established integral representations of these solutions, best known as the 'path integral' 
representation. Explicit results were obtained for (p 3 , ip 4 , if 3 + ip 4 and the charged scalar 
field theories. Subsequently, we studied the 'renormalization' of such theories in zero di- 
mensions, which is equivalent to counting diagrams with restrictions imposed on the type 
of diagrams considered, for instance diagrams without any tadpoles, self-energy insertions 
or vertex insertions. We were able to get explicit results for the dependence of the bare 
quantities such as the mass, the coupling, and the tadpole counter terms, on the renormal- 
ized (physical) coupling constant. Examples of interesting observations are the facts that 
in ip 4 theory, the bare mass exhibits a zero at a finite value of the renormalized coupling 
constant g = g+ (Eq.fl66D), whereas at g — > 2 — e the bare coupling becomes strong and 
the mass squared becomes large and negative. Yet in both cases the 'physical' connected 
Green's functions remain finite and calculable. 
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Appendix 

Consider general ^-theory, and suppose that all but one renormalization conditions have 
been implemented through functions z k (g, z p ), k = 1, 2, ..,p — 1 of two variables z p and g, 
like in Eq. ([53"l) and Eq.(^TJ). This means that we are considering a theory with an action 

S(g,z p] (p) = —(p p + }^— — up . 

P ' k=l 

Let T be a contour in the complex yj-plane, such that Re(p n — > oo at the endpoints, and 
define 

Z n (g,z p ):= J d(f(p n exp{-S(g,z p ; ip)} . 

Such an integral is not defined for all complex values of gz p . Let gz p = \gz p \e ipri and denote 
by e~ iv T the contour that is obtained from F by clockwise rotation over rj. For complex 
values of gz p , we define 



Z„{g,z p ):= / dupuo n exp{-S(g,z p ;uj)} 

Je- ir T 



k=l ^ 

Integrals of this type can easily be calculated to high precision by numerical integration. 
One just has to choose T such that it goes through one or more saddle points, so that the 
integrand oscillates as little as possible. 

To formulate the renormalization problem further, let us denote the connected moments 
by C n , so 

n — 1 n — ^ 2 r 12 n _ ^ 3 in n n?> 

•--l — TT i ( - J 2 — °1 ) °3 — ~ 00 2 0i — O x , . . . 

and so on. The problem is to solve z p as function of g from the implicit function equation 

C P (g, z p ) = -g , 

which represents the final renormalization condition. This equation can be solved numeri- 
cally. Given a value of g, we have to find the zero of the function 

F(zp) := Cp(g,z p ) +g , 

which can be found using Newton-Raphson iteration 

z ^z F{Zp) 
P P F'(zp) ■ 



36 



By making small steps in the value of g, the solution z p (g) on a curve in the complex g- 
plane can be determined. At the start of each iteration, the question arises of which initial 
value of z p to choose, and the obvious answer is to choose the final value of the previous 
iteration, which should lie close the the new final value if the steps in g are not to large. 

As a check one can look whether the results obtained with this method satisfy (numer- 
ically) the available differential equations for z p (g) (Eq. (|55|) and Eq. (|63|) ). 
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